Intermediate liquid-solid regime near the quantum melting 
of a two dimensional Wigner molecule 
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For intermediate values of the Coulomb energy to Fermi energy ratio r s , the ground state of a 
few spinless fermions confined on a two dimensional torus is the quantum superposition of a floppy 
Wigner molecule with delocalized vacancies and of a Fermi liquid of intersticial particles. This raises 
the question of the existence of an unnoticed liquid-solid phase between the Fermi liquid and the 
Wigner crystal for fermionic systems in two dimensions. 

PACS: 73, 73.20.-r, 73.21.-b, 67.80.-s 
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An outstanding question in quantum condensed mat- 
ter is to know if the solid and the fluid can coexist in an 
intermediate phase separating the solid from the liquid. 
Such a phase was suggested if the zero point motions 
of certain defects become sufficient to form waves propa- 
gating inside the solid. The defects of the crystalline solid 
can be simple vacancy-intersticial pairs or more complex 
excitations. Their statistics coincide with the statistics of 
the particles out of which the solid is made. For bosons, 
they may form a condensate, giving rise to a superfluid 
coexisting with the solid. This supersolid phase is dis- 
cussed in certain bosonic models [pj. For fermions, the 
defects may form a Fermi liquid B coexisting with the 
solid, such that the system is neither a solid, nor a liquid. 
Two kinds of motion are possible in it; one possesses the 
properties of motion in an elastic solid, the second pos- 
sesses the properties of motion in a liquid. 

We study such a possibility for electrons created in a 
two dimensional (2d) heterostructure. Our motivation is 
threefold. First, it becomes possible to create very dilute 
2d gases of electrons or holes in field effect devices and 
to vary by a gate the carrier density n s such that the 
Coulomb energy to Fermi energy ratio r s cx n s ~ x l 2 can 
reach high values where the charges become strongly cor- 
related. The studied devices use doped semi-conductors 
Pj (Si-Mosfet, Ga-As heterostructures, Si-Ge quantum 
wells) or undoped organic crystals |^| . Since charge crys- 
tallization is expected for r s « 37 in a clean system 
and at lower r s in the presence of impurities [Q,^), to 
study in those devices how one goes from a Fermi liquid 
(low r s ) towards a Wigner crystal (large r s ) becomes pos- 
sible. Second, the observation Q of a metallic behavior 
at low temperatures in those systems raises the question 
of a possible intermediate phase, which should be neither 
a Fermi system of localized particles (Anderson insula- 
tor), nor a correlated and rigid solid of charges (pinned 
insulating Wigner crystal). Third, an unexplained inter- 
mediate regime was numerically observed in small 
disordered clusters. Those studies give a first threshold 
rf above which the Hartree-Fock (HF) approximation 
does not describe [ fl0| the persistent currents driven by 
an Aharonov-Bohm flux <b. For r s < , the local cur- 



rents are oriented at random by the impurities, while 
for r s > rp they align along the shortest direc- 

tion enclosing <j>. 



is followed by a second threshold 
where charge crystallization occurs and the persistent 
currents vanish. This intermediate regime characterizes 
the ground state and the low energy excitations and dis- 
appears when the excitation energy exceeds |j| the Fermi 
energy. Both numerics and experiments give an unex- 
plained behavior restricted to temperatures smaller than 
the Fermi temperature for similar intermediate values of 
the ratio r s . 

Since the low energy levels do not obey || Wigner- 
Dyson statistics, the existence of non chaotic low energy 
collective excitations due to the interplay between the 
kinetic energy and the Coulomb repulsion in the clean 
limit can be suspected. For this reason, we have stud- 
ied the same cluster than in Refs. |^|| without random 
substrate and we have observed for intermediate values 
of r s a liquid-solid regime of a type conjectured by An- 
dreev and Lifshitz. The transfer of a charge from a crys- 
tal site towards some intersticial site creates a vacancy- 
intersticial pair. This costs an energy eo, but the pair can 
tunnel to give rise to a band of delocalized excitations of 
width to- Near the melting point of the crystal, to may 
exceed eo and the Wigner solid may co-exist with a liq- 
uid of delocalized defects. This picture is supported by 
an analysis of the ground state (GS) using a combination 
of Slater determinants (SDs) built out from plane waves 
and from site orbitals. The plane wave SDs are given 
by some excitations of the non interacting system and 
correspond to an excited liquid. The site SDs describe 
the Wigner solid molecule which dominates at large r s . 
For rf 9.31 < r s < rf » 27.93 in the studied cluster, 
the GS is given by a quantum superposition of an excited 
liquid and of a Wigner solid. The structure of the low en- 
ergy excitations, the errors made assuming certain trun- 
cations of the Hilbert space and the responses to small 
perturbations confirm the existence of two thresholds rf 
and rf between which unusual behaviors occur. 

We consider N spinless fermions free to move on a 
square L x L lattice with periodic boundary conditions 
(BCs). As in Refs. gg, we take N = 4 and L = 6, the 
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size Nh = 58905 of the Hilbert space being small enough 
to exactly diagonalize the Hamiltonian: 



H = -tJ2 4c j + U^ nin i 



<«,j> 



(1) 



using Lanczos algorithm, c] (a) creates (destroys) a spin- 
less fermion in the site i, t is the hopping term between 
nearest neighbors and ry is the shortest inter-particle 
distance in a lattice with periodic BCs. rii — c[ci and 
the interaction strength U yields a Coulomb energy to 
Fermi energy ratio r s = U / (2ty/im s ) for a filling factor 
n s = N/L 2 = 1/9. 

The system remains invariant under rotation of an- 
gle 7r/2 and under translations and reflections along the 
longitudinal x and transverse y directions. Invariancc 
under translations implies that the momentum K is a 
good quantum number which remains unchanged when 
U varies. The symmetries imply that the states are four- 
fold degenerate if K ^ and can be non degenerate if 
K = 0. 

When U = 0, the states are Nh plane wave SDs 
4(4)4(3)4(2)4(1)1° > (4( P ) creating a particle in a 
state of momentum k(p) — 2ir(p Xl p y ) / L), |0 > being 
the vacuum state and p XiV = 1, . . . ,L). Ki and Di being 
the momentum and the degeneracy of the states of en- 
ergy Ek, the GSs (E (U = 0) = -13*, #0 ^ 0,D = 4) 
are followed by two sets of excitations {E\{U = 0) = 
-12t,A = 25) and (E 2 (U = 0) = -11*, D 2 = 64). 
We denote \K {/3) > (/3 = 1,...,4) the SDs of energy 
— 13* and of momenta K = (0,±7r/3) and (±7r/3, 0) and 
\Ki((3) > the 4 SDs of energy —12* and of momentum 
K\ = 0. For the non interacting system, the \Kq((3) > 
are 4 orthonormal GSs and the K\ (/?) are 4 orthonormal 
first excitations, corresponding to a particle at an energy 
—4* with fc(l) = (0,0), two particles at an energy —3* 
and a fourth particle of energy —2t with momenta such 
that £j= 2 fc (j) = 0. One has fc(2) = (0,±tt/3), fc(3) = 
(±tt/3,0) and fc(4) = (t^A or fc(2) = (0,=Ftt/3), 

fc(3) = (±tt/3,0) and fc(4) = (=prA ±tt/3). 

When * = 0, the states are Nh Slater determinants 
cjc]c|.c||0 > built out from the site orbitals. The con- 
figurations ijkl correspond to the Nh different patterns 
characterizing 4 different sites of the L x L square lat- 
tice. If we order the configurations by the smallest dis- 
tance d between two sites, Nd denoting the number of 
configurations with inter-particle spacings larger than d, 
one has Nj. = 27225, N^ = 9837, N 2 = 2709, N^ = 81 
configurations out of 58905 configurations. The N^ con- 
figurations have the smallest electrostatic energies, and 
contain 9 square configurations \Dj > (I = 1, ... ,9) of 
side a = 3 (energy E (t = 0) « 1.8017), 36 parallelograms 
of sides (3, vlO) (energy w 1.85J7) and 36 other parallel- 
ograms of sides (a/10, \/l0) (energy w 1.97J7). Ordering 
the site configurations by increasing electrostatic energy, 
those 81 configurations are followed by 144 configurations 



obtained by moving a single site of a square configura- 
tion by one lattice spacing (deformed squares of energy 
S3 217). 




FIG. 1. As a function of r 3 , low energy part of the spectrum 
exhibiting a level crossing at rf . Inset: two first level spac- 
ings Ai/t (dashed) and A2/* (dotted) and the perturbative 
result 3DrJ 3 (dot-dashed). 



The low energy part of the spectrum is shown in Fig. 
[l] as a function of r s . If we follow the 4 GSs Eq{t s = 0) 
(Kq 0), one can see a first level crossing with a non 
degenerate state (iio = 0) which becomes the GS above 
rf , followed by two other crossings with two other sets 
of 4 states with Kj 7^ 0. When r s is large, 9 states 
coming from E\(r s = 0) have a smaller energy than 
the 4 states coming from E (r s = 0). The degenera- 
cies ordered by increasing energy become (1,4,4,4) in- 
stead of (4, 25, 64) for r s = 0. Since the degeneracies are 
(9,36,36) when * = 0, these 9 states give the 9 square 
molecules |Dj > when r s — > 00. The centers of mass Ri 
of the \Oj > are located on a periodic 3x3 square lattice. 
When r s is large, one has a single massive molecule free 
to move on this restricted lattice, with a hopping term 
T oc trj 3 and quantized Ki(I) — 2irpi/3 longitudinal 
and K t (I) = 27rp t /3 transverse momenta (pi t t — 1,2,3). 
This gives 9 states of different kinetic energies given by 
-2T{cosKi(I) + cosK t (I)). The kinetic part of the low 
energy spectrum is then — AT, —T, +2T with degeneracies 
1, 4, 4 respectively. This structure with two equal energy 
spacings Ai and A2 appears (inset of Fig. ^ when r, 
is larger than the crystallization threshold ■ Above 
, to create a defect in the rigid molecule costs a high 
energy available in the lQ th excitation only. 

To describe large r s , one can use degenerate perturba- 
tion theory and study how the degeneracy of the 9 |dj > 
is removed by terms oc U/t oc r^ 1 . One obtains for the 9 
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low energy states: 

^ = E D - 2^(cos K t (I) + cos tf t (I)). (2) 

£d = Ar s + B/r s + C/rl and T/i = D/rj (A w 2.13, 
5 w -70.81, C w -18763.48 and L> w 3463.97). 
i?_D comes from the small elastic vibrations of the rigid 
molecule while 8T is the band width of its zero point 
fluctuations. 



Below rf, each of the 4 GSs \V%(r a ) > with K ^ has 
still a large projection 




FIG. 2. Errors AE {N d ) = (E {N d ) - E ) /t as a function 
of r s : thick line (d — 1), dotted line (d — V2), circle (d = 2), 
square (d = v5). 

The two thresholds rf" and can be also detected 
if one calculates the GS energy Eo(Nd) of the Hamilto- 
nian truncated onto the Nd site SDs having a minimum 
inter-particle spacing > d. The difference AEo(Nd) = 
(Eo(Nd) — Eo)/t between the GS of the truncated Hamil- 
tonian and the exact GS is given in Fig. ^. AEq(Ni) be- 
comes negligible above rf where the probability to have 
two nearest neighbor particles becomes weak. AEq (AOg) 
coincides with AE (n) = E (t = 0) - E . The curve 
AE (N 2 ) is very remarkable. Above rf , AE (N 2 ) 
is identical to AEq(0), while for rf < r s < r^ , 
AEq(N2) « 3 independently of r s . This suggests that the 
GS for intermediate r s is composed of an elastic molecule, 
which can be projected onto the N 2 site SDs adapted to 
describe it, plus "excitations" having a nearly interac- 
tion independent kinetic energy — 3i, i.e. the energy 
of a particle at the Fermi surface of the non interacting 
system. 

To understand further the nature of the intermediate 
GS, we have projected the GS wave functions l^oC^s) > 
over the two eigenbases valid for U/t = (Fig. |^ upper 
left) and for t/U = (Fig. [| upper right) respectively. 



Po(r s ) = J2\ <xS, o(r s )\KoW> 

13=1 



(3) 



over the 4 non interacting GSs. There is no projection 
over the 25 first excitations and a small one P2{r s ) over 
the 64 second excitations of the non interacting system. 
Above rf , the non degenerate GS with Kq = has a 
large projection 
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Ei 
0=1 



< * (r.)|ifi(/3) > 



(4) 



which is equally distributed over the 4 excitations 
\K\{0) > of momentum K\ = 0. Its projections onto 
the 4 \Kq((3) >, the 21 other first excitations and 64 sec- 
ond excitations of the non interacting system are zero. 
One concludes that a large part of the system remains 
an excited liquid above rf , described by 4 equal projec- 
tions onto the 4 \K\(/3) >. Those projections decrease as 
r s increases. 

The modulus A(U) = \ < a\^ (r s ) > \ of the ampli- 
tude of the projection of the GS |\l/o( r s) > over a sin- 
gle square configuration |D > is given in Fig. || (upper 
right). A(D) does not depend on the chosen square and 
reaches 1/3 for large r s , i.e. the value where the total 
projection over the 9 squares is equal to one. A(0) lin- 
early increases when rf < r s < r^ and gives r s ps 37 
(the accepted value for 2d Wigner crystallization) by ex- 
trapolation. One concludes that a small but increasing 
part of the system begins to be a Wigner square molecule 
at rf , described by 9 equal projections onto the \Dj >. 
The site SDs and plane wave SDs are not orthonormal. 
After re-orthonormalization, the total projection p t of 
|4"o(^s) > over the subspace spanned by the 4 \Ki(f3) > 
and 9 |Dj > and P t over the subspace spanned by the 4 
\K\{0) > and 225 site SDs of lower electrostatic energies 
(9 squares, 36 + 36 parallelograms, 144 deformed squares) 
are given in an inset of Fig. |^ (upper right). One can see 
that l^o^s) > mainly remains inside this small part of 
the large Hilbert space for intermediate r s . 

We consider now the GS response to small perturba- 
tions. The first one consists in piercing the 2d torus by 
an infinitesimal positive flux <f> (periodic tranverse BCs, 
t — ► iexp(i(/)/L) for longitudinal hopping only, (j) = tt 
corresponding to anti-periodic longitudinal BCs). The 
coefficients a(r s ) and b(r s ) (Kohn curvature) of the ex- 
pansion E (r 8 ,(j>) w E Q (r s ,0) + a(r s )<f> + b(r s )(j) 2 /2 are 
given in Fig. ^ (lower right). When r s = 0, (f> removes the 
fourfold degeneracy of Eo, a = -y/Et/6 and b = 7t/36. 
When r s is large, the substitution Ki(I) — » Ki(I) + 2<p/3 
in Eq. || gives a = and b w 8DtrJ 3 /9. An infinites- 
imal positive flux <j> gives rise to a persistent current 
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Ii = —dE /d(j) — —a when r s < rf while the GS curva- 
ture b exhibits a smooth crossover between two regimes 
around (lower left inset of Fig. ft) . 
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FIG. 3. As a function of r a : Upper left: GS Projections 
Po (dotted), Pi (thick) and P2 (dot-dashed) onto plane wave 
SDs; Upper right: GS Amplitude A(D) (thick) over a single 
square site SD (the dot-dashed line reaches 1/3 at r s m 37); 
Inset: GS Projection p t (dashed) and Pt (dotted) over com- 
bined plane wave and site SD re-othonormalized bases; Lower 
left: a(r„)/t (dotted) and b(r s )/t (thick) characterizing the 
GS response to an infinitesimal flux; inset : b(r 3 )/t (dashed) 
and SD/9r~ [i (dotted); Lower right: GS density D p (r s ) at the 
pinning site with V p /t — —0.01 (dotted), —0.05 (dashed) and 
-0.1 (dot-dashed). 

The second perturbation consists in introducing a weak 
pinning well (negative potential V p ) at a single lattice site 
p. The GS density 



D p (r s ) =< *o(r»)l4cpl*o(r. 



> 



(5) 



at the site p is shown in Fig. || (lower right). If V p — 0, 
D p (r s = 0) = 1/9. A weak negative value of V p yields a 
larger value for D p (r s = 0). When one turns on the inter- 
action, the effect is first enhanced and enhanced Friedel 
oscillations around p can be observed. At rf, D p drops 
and the interacting GS has a weaker response to a weak 
pinning well than the non interacting GS, as shown by 
the curves D p (r s ) and by the weaker Friedel oscillations 
observed around p. When r s is large, D p increases again 
and the Wigner molecule is pinned. This surprizingly 
weak response for intermediate r s suggests that the sys- 
tem may very weakly respond to the presence of weak 
impurities. 

From the study of the GS projections emerges the con- 
clusion that a minimal description of the intermediate GS 
requires to combine the two limiting eigenbases. Instead 



of taking separately many low energy excitations of the 
U/t = limit or of the t/U — limit, one can use the 
subspace spanned by the 9 ]□/ > and the 4 l-ft'i(O) > 
for having a large fraction of | 'Jo (?";;) > for intermediate 
r s . In this sense, the GS is neither solid, nor liquid, but 
rather the quantum superposition of those two states of 
matter, as conjectured in Ref. EL This suggests possi- 
ble improvements of the trial GS to use for intermedi- 
ate r s in variational quantum Monte Carlo approaches 
Instead of using Jastrow wave functions improving 
the plane wave SDs for the liquid or the site SDs for 
the solid, it will be interesting to study if a combina- 
tion of the two describing a solid-liquid regime will not 
be better for intermediate r s . This will confirm that an 
unnnoticed intermediate solid-liquid phase does exist in 
the thermodynamic limit for fermionic systems in two 
dimensions. This might help to explain the 2d metal ob- 
served in field effect devices. Eventually, our results for 
spinless fermions might also give a natural explanation 
for the competition between the Stoner ferromagnetism 
and the Wigner antiferromagnetism observed (l2| for in- 
termediate r s when the spin degrees of freedom and the 
disorder are included. 
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EU. 
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